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Abstract 

I We study the production of the stable six-quark H dibaryon via the (K~,K + ) reaction on a 12 C target 

O ; within a covariant effective Lagrangian model. We follow a factorization approximation, in which the full 

SZ production amplitude is written as a product of the amplitudes for the K + p —* K + + S and H + p — > H 

i \ processes. The K + "BT production vertex is described by excitation, propagation and decay of A and X 

ly-j " resonance states in the initial collision of a K meson with a target proton in the incident channel. The 

t>" 

t^. ■ parameters of the resonance vertices are taken to be the same as those determined previously by describing 

o: 

the available data on total and differential cross sections for the p(K , K + )E reaction within a similar 
I model. The 3 + p — > H fusion process is treated within a quark model where the H dibaryon is considered 

t-H ■ 

as a stable particle. For the K + meson angle fixed at 0° , the H production cross-section is found to be about 
2.7 pb/sr for H mass just below the AA threshold at a K beam momentum of 1.67 GeV/c. This is an order 
of magnitude larger than the value for this quantity reported earlier in calculations performed on a He target 
using a different model for the cascade hyperon production. We have also calculated the beam momentum 
dependence of the H production cross section and the energy spectrum of the emitted K + meson. 

PACS numbers: 13.75.Jz, 14.20.pt, 25. 80.Nv 
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Within the quark-bag model, the H dibaryon, a six-quark [two up (u), two down (d) and two 
strange (s)] state with spin-parity J" = + , and isospin 1 = 0, was predicted to be a stable system 
with a mass about 80 MeV below the AA threshold, some 35 years ago yj]. Later calculations that 
included the center-of-mass (CM) [2] and pionic-cloud |3|] corrections within this model, predicted 
this state to be much less bound or even unbound. Since then a lot of experimental effort has gone 
into searching for the H dibaryon (see, e.g., the review Q4£| for references up to the year 2000 and 
Refs. [5] and [6] for more recent investigations done by Japan's National Laboratory for High 
Energy Physics (KEK) and STAR at RHIC collaborations, respectively). These studies have led to 
the conclusion that the existence of this system as a deeply bound object is highly unlikely. At the 
same time, the observation of the double-A hypernuclues AA He (NAGARA event) and the precise 
determination of its binding energy at KEK in the experiment E373 [7] has put a lower limit of 
2.2237 GeV to the //-dibaryon mass at a 90% confidence level, which is just about 7.7 MeV below 
the AA threshold. 

The interest in the //-dibaryon has been revived by the recent lattice quantum chromodynamics 
(LQCD) calculations of different groups. The NPLQCD [8|] and HAL QCD |9J] collaborations 
have reported that the H particle is indeed bound at somewhat larger than physical pion masses. 
However, extrapolations of the calculations of these groups to the physical pion mass region sug- 
gest J^ lll| that this particle is likely to be either very loosely bound or an unbound state near the 
AA threshold. In a very recent chiral constituent quark model calculation [jl2|], the value extracted 
for the binding energy of the H particle has been found to be compatible with the restrictions im- 
posed by the NAGARA event. These results together with the previous experiments [5] that give 
a hint on the existence of a possible //-dibaryon resonance, have led to a proposal to look for this 
particle in a future experiment IU3I1 at the Japan Proton Accelerator Research Complex (JPARC) 
using a high intensity K~ beam. This is expected to answer the long standing question about the 
existence of the //-dibaryon. Furthermore, with the measurement of the exclusive 'Br production 
in the yp —* K + K + E~ reaction at the Jefferson Laboratory yjj], a possibility has been opened for 
producing the //-dibaryon with the photon beam. 

The (K~,K + ) reaction leads to the transfer of two units of both charge and strangeness to the 
target nucleus. Thus this reaction is one of the most promising ways of studying the production of 
S = -2 systems such as S hypernuclei and the //-dibaryon. Recently, the production of cascade 
hypernuclei via the (K~,K + ) reaction on nuclear targets, has been investigated within an effective 



Lagrangian model fl 1 3M - This is a new approach, where the K + E production vertex is described by 



excitation, propagation and decay of A and 2 resonance intermediate states in the initial collision 
of the K meson with a target proton in the incident channel. The S hyperon gets captured into 
one of the nuclear orbits leading to the formation of the E~ hypernucleus. The parameters at the 
resonance vertices were determined by describing the available data on total and differential cross 
sections for the elementary process p(K~ , K + )E~ within a similar approach Jl^L where contribu- 
tions were included from the ^-channel and w-channel diagrams, which have as intermediate states 
A and 2 hyperons together with eight of their three-and four-star resonances [A(1405), A(1520), 
A(1670), A(1810), A(1890), 2(1385), 2(1670) and 2(1750)] with masses up to 2 GeV. It was ob- 
served in this study that the total cross section of the p(K~ , K + )E~ reaction is dominated by the 
contributions from the A(1520) (with L u = D 0i ) resonance intermediate state. The region for 
beam momentum (pk) below 2.0 GeV/c was shown to be dominated by contributions from the 
5-channel graphs - the w-channel terms are dominant only in the region p K - 1 2.5 GeV. 

In calculations of the S~ hypernuclear production cross sections, one requires additionally the 
bound state spinors for the proton hole and E~ particle bound states. These were obtained by solv- 
ing the Dirac equation with vector and scalar potential fields having Wood-Saxon shapes. Their 
depths were fitted to the binding energies of the respective states. In Ref. Ill5ll bound state spinors 



obtained in the quark-meson coupling (QMC) model II 1711 were also used. The cross sections for 



the hypernuclear production were found to be quite different from those calculated previously in 



Ref. 



18] 



In this paper, we describe the production of the //-dibaryon by the (K~,K + ) reaction on a 
12 C target using a similar approach. The basic production mechanism considered in our work is 
depicted in Fig. 1, where the reaction proceeds in two steps. In the first step the K + E production 
takes place by following the process as described above while in the second step the S~ hyperon 
fuses with another proton of the residual nucleus to form the H dibaryon. A similar method was 



also used earlier in Refs. Bfl in calculations of the production of this particle via (K , K + ). 
However, there are two important differences between those calculations and the present work, 
(i) We use the proper three-body phase space and kinematics for the K~+ n C — » K + + //+ 10 Be 
process, (ii) We calculate the K + E production amplitude by employing the method described in 
the previous paragraph that was applied in our calculations of the S~ hypernuclei [15], while these 
authors used a parameterization of the sparsely available experimental zero degree differential 
cross-section for the p(K~ , K + )E~ reaction and multiplied them by the target wave function. The 
earlier calculations were limited to a very light 3 He target nucleus, whereas we apply our method 
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FIG. 1. (color online) Graphical representation of the model used to describe the p{K~ , K + )E~ (Fig. la) 
and A(K ,K + )HB reactions (Fig. lb), where A represents the target nucleus while B = (A-2) the residual 
nucleus. In Fig. la, A* and E* represent the A and £ resonance states, respectively. 

to 12 C. This makes it possible to compare our cross sections directly to the existing experimental 
results and to make predictions for future measurements. 

It should be remarked that there may be other processes also through which the //-dibaryon 
production can proceed. These are termed higher order processes in Ref. [20], where it was shown 
that their contributions are not expected to be large. Like these authors, we have neglected such 
diagrams in our study. 

The general formula for the invariant cross section of the K+A^>K + +H + B reaction is 



written as(see,e.g., Ref. Il22ll ). 
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where A f t represents the total amplitude, P, and P f the sum of all the momenta in the initial and 
final states, respectively, and m H , m A and m B the masses of //-dibaryon, and nuclei A and B, 
respectively. The cross sections in the laboratory or CM systems can be written from this equation 
by imposing the relevant conditions. It may be noted that summations over final spin states and 
average over initial spin spates are implied in Ay>. 

Following the factorization approximation of Ref. [20J], the total amplitude A j, is written as the 
product of the amplitudes for the processes K~ + p\ — » K + + ET (M Pk _^ Pi ^_), and S~ + p 2 —* H 
(F P2 ,p 3 -) (see Fig. lb). We write 



/d A p\ C d 4 p2 
(2nf J llnf 5 ^ 1 + Pk ~ ~ Py '^ Py ' ~ Pk+ ~ Pz-WP" ~ P*- ~ P^ 

x [Tj M pk-,pup s -( K ~ +Pi ~* K + +S-)]F P2 ,p s -(S-+p 2 ^ H), (2) 
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where p K -, Pk + , Py*, and p s - are the four momenta of the incoming and outgoing kaons, inter- 
mediate resonance state and the E~ hyperon, respectively. The amplitude M PK _ tPuP7 ._, where the 
summation is done over all the resonance intermediate states Y* as described above, has been 
determined by following the method discussed in Ref. Ill 511 . The effective Lagrangians, the cor- 
responding coupling constants and the form factors for the resonance-kaon-baryon vertices, the 
propagators for the intermediate resonances and the bound state and free-space wave functions for 
the bound proton p\ and the intermediate S~ hyperon, respectively, as used in the calculations if 
the amplitude M, are discussed in the following. 

The effective Lagrangians for the resonance-kaon-baryon vertices for spin-^ and spin-| reso- 
nances are taken as 



M 



£ 



gKBR 



KBR 3/2 



3/2 



m K 



¥ Rv d^K^B + h. c, 



(3) 
(4) 



with M = {m R + m B ), where the upper sign corresponds to an even-parity and the lower sign to 
an odd-parity resonance with B representing either a nucleon or a S hyperon. The operator Y is y 5 
(1) for an even- (odd-) parity resonance. The parameter ;f controls the admixture of pseudoscalar 
and pseudovector components. The value of this parameter is taken to be 0.5 for the A* and Z* 
states, but zero for A and S states, implying pure pseudovector couplings for the corresponding 
vertices, in agreement with Refs. 11221. 12311 . It may be noted that the Lagrangian for spin-|, as given 
by Eq. (4), corresponds to that of a pure Rarita-S winger form which has been used in all previ- 
ous calculations of the hypernuclear production reactions within a similar effective Lagrangian 
mod el Il23l - l25|l . The values of the vertex parameters were taken to be the same as those given in 
Ref. Ua. 



Similar to Refs. [15, 



1611 . we have used the following form factor at various vertices, 

,! 4 



F m (s) = 



(5) 



/l 4 + 0-m 2 ) 2 ' 

where m is the mass of the propagating particle and A the cutoff parameter, which is taken to be 



1.2 GeV, which is the same as that used in Refs. [115 



w. 



The two interaction vertices of Fig.Q]are connected by a resonance propagator. For the spin- 1/2 
and spin-3/2 resonances, the propagators are given by 

KjnP 1 + m Rl/z ) 



Vr 



1/2 



P 2 -(m Rl/2 -ir Rl/2 /2) 2 ' 



(6) 



1/2' 



and 



i( nP +m Rm ) 



p 2 -(m Ri/2 -iY Ry2 /2y 



respectively. In Eq. © we have defined 

P »y = r _ - t^p"p v + t~ — (p"/ - pV) ■ ( 8 ) 

3 3m^ a/2 3m R3/2 

In Eqs. © and ©, Y Rl/2 and r Ra/2 define the total widths of the corresponding resonances. We 
have ignored any medium modification of the resonance widths while calculating the amplitude 
M, as information about them is scarce and uncertain. 

The bound proton wave function [i]j(p\ )] is a four component Dirac spinor that is the solution of 
the Dirac Equation for a bound state problem in the presence of external scalar and vector potential 
fields. This is written as 



( mw;i /2j (pi) 

-ig(h)y™ J l/2j (Pi), 



(9) 



1/2/ 

where f(h) and g(k { ) are the radial parts of the upper and lower components of the spinor ip(pi). 
In this equation, J/"j J /2 y are the coupled spherical harmonics 

y'n/ij = < im { \l2pi\jmj > Y tm {p\)X^ ( 10 ) 

and £' = 2£ - j with I and j being the orbital and total angular momenta, respectively. Y represents 
the spherical harmonics, and^ the spin space wave function of a spin-^ particle. 

We assume that the nucleon bound state has a pure single-hole configuration with the core 
remaining inert. The spinors in momentum space are obtained by Fourier transformation of the 
corresponding coordinate space spinors, which are the solutions of the Dirac equation with po- 
tential fields consisting of an attractive scalar part (V s ) and a repulsive vector part (V v ) having 
a Woods-Saxon form. For fixed geometry parameters (radius and diffuseness) we search for the 
depths of these potentials to reproduce the binding energy of the respective state. For the p\ proton 
(see Fig. lb) this state has the quantum numbers Ip^j with a binding energy of 15.96 MeV. The 
searched depths were 382.6 MeV and -472.3 MeV, respectively for the fields V v and V s with the 
radius and diffuseness parameters of 0.983 fm and 0.606 fm, respectively for both. It should be 



added here that these spinors were the same as those used in Refs. [25] the production of the A 
hypernuclei via the (y, K + ) reaction on a 12 C target. 
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FIG. 2. (color online) Differential cross section do-jdQ.^ at 9 K + - 0° for H production in the 
l2 C(K~ , K + )H l0 Be reaction at the beam momentum of 1.67 GeV/c, as a function for //-dibaryon mass. 
The AA and 1.N thresholds are shown by inverted arrows as indicated. 



The free-space spinors for the 3 hyperon are written as 



y(ps-) = 8(p s - -E s -) 



Es- + tn~.- 



<r-Pn- 



(ID 



After having established the vertices and the corresponding coupling constants as well as the 
forms of the propagators, the amplitudes M can be written by following the well known Feynman 
rules and can be numerically evaluated using the bound state and the free-space spinors. We have 
used a plane wave approximation to describe the relative motions of the kaons in the incoming 
and outgoing channels. However, the distortion effects are partially accounted for by introducing 
reduction factors to the cross sections as described in Ref. tl5|]. It may be noted that these factors 
correspond to absorption effects on only the K~ and K + wave functions. There could still be the 
distortion effect on the H - B relative motion in the final channel which is ignored here. Since 
all of our calculations are carried out in momentum space, they include all the nonlocalities in the 
production amplitudes that arise from the resonance propagators. 



For calculating the amplitude F p ^ p __, we follow the same procedure as described in Ref. Q20Q . 
In this method, the //-dibaryon is treated as a bound particle with a mass m H , however its six quark 
structure is taken into account. This implies that three-quark internal structure of the hyperon 
and the proton p 2 (see Fig. lb) also have to be invoked as the formation of H is thought of in terms 
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of the fusion of two three-quark bags (S - and p 2 ). The proton p 2 is supposed to be picked up 
from the same orbit as p\ and has a binding energy of 1 1.22 MeV. The amplitude F is calculated 
by taking the overlap of the internal wave functions of H, E~ and p 2 , which are described by a 
Gaussian approximation (see, Ref. Jl^l ). The final result for the amplitude F is given by 

w t~±. m r ( 2R p Rh \ 3 ( 2R ^ Rh \ 3 ( R h\ 3 ' 4 

f Pp2 , Ps _(* +P2 ^H) = ^(-r^) (^-) ( 3^) 

R 2 

Xexp[--^(p 2 -p s -) z ], (12) 



where the factor r arises from the color- flavor- spin recoupling as defined in Ref. 112011 . Its value 
is VI /20. The values of the oscillator parameters, R p , and R H have been taken to be 0.83 fm, 
0.73 fm and 0.95 fm, respectively. The chosen value of R p reproduces the root mean square (rms) 
radius of the proton, while the value of R~. comes from the quark-bag model relation between 
the proton bag radius and that of the S - . The bag radius of the //-dibaryon is about 20% larger 
than that of the proton (see, eg. Refs. iQjIs]). Therefore, we have increased the value of R H over 
R p accordingly. We do not use the approximation R p = i? s = R H as done in Refs. 
should be remarked here that in deriving Eq. [T2l the normalizations of the baryon wave functions 
are made consistent with those used for them in the amplitude M. 

The method discussed above has been used to calculate the cross-sections for the 12 C(K~ , K + )H w Be 
reaction. In Fig. 2, we show the results for the differential cross section dcr/dQ. K + at the K + angle 
of 0° and the beam momentum (pk ) of 1.67 GeV/c as a function of the rest mass of //-dibaryon. 
We see that the cross sections decrease uniformly as m H approaches the AA threshold (indicated 
by an arrow) where we stopped the calculations because our method treats H as a bound particle. 
It should be noted that near this threshold our cross section is about an order of magnitude larger 

n El 

than that reported in Ref. H20I1 at the beam momentum of 1.8 GeV/c on a He target. In Ref. [5J, 
the production cross section of H with a mass range between AA and 2/V thresholds (also indi- 
cated in Fig. 2 by an arrow) has been found to be 2.1 + 0.6 (stat.) +0.1 (syst.) pb/sr at a 90% 
confidence level, in a measurement of the l2 C(K~ , K + )AAX reaction at the beam momentum of 
1.67 GeV/c where K + meson was confined mostly in the forward directions. We see in Fig. 2 that 
our cross-section at the AA threshold is comparable to this value. 

In Fig. 3, we show the beam momentum dependence of the dcr/dQ,K + of the same reaction as 
that of Fig. 2. We note that the cross section peaks near the p K - value of about 1 .3 GeV/c, which is 
approximately 0.50 GeV/c away from the production threshold for this reaction. This is similar to 
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p K - (GeV/c) 

FIG. 3. (color online) Differential cross section dcr/d£lK + at 9k+ = 0° for the same reactions as in Fig. 2 
reaction as a function of the beam momentum px- corresponding to a Mh of 2.224 GeV. 



the situation in the case of the zero degree differential cross section of the elementary production 
reaction p(K~ , K + )E _ , where the peak occurs near a value of px-, which is away by a similar 



amount from the corresponding production threshold (see Fig. 4 of Ref. [15]). On the other hand, 
the distribution in Fig. 3 is broader as compared to that of the elementary reaction. This can be 
attributed to the three-body phase-space factor (see the discussions below) and also partly to the 
Fermi motion of the protons in the target nucleus. 

If one compares the magnitude of the peak cross section in Fig. 3 with that of the elementary 
production reaction, one notices that the H production cross section is smaller than that of the 
p(K~ , K + )E" reaction by at least an order of magnitude. This reflects the fact that H production 
occurs only in a limited region of S~ - p relative momenta. This figure further shows that in 



contrast to conclusions of Refs. H2011 and H21H . using a K~ beam of lower momentum, around 1.4 
GeV/c, provides the maximum probability of producing //-dibaryon via the (K~,K + ) reaction on 
a 12 C target. It is expected that use of a heavier target like Cu (which is supposed to be employed 



in the experiment to be conducted in the proposal of Ref. 111311 ) will push the peak in the cross 
section to even lower beam momenta because the corresponding production threshold would be 
lower than that on a 12 C target (the threshold beam momentum for the 63 Cu target is 0.701 GeV as 
compared to 0.735 GeV for the 12 C corresponding to a M H of 2.224 GeV). 
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FIG. 4. (color online) (a) Double differential cross section dcr/dD.K+dpK+ at 8k+ = 0° for the 
l2 C(K~, K + )H m Be reaction as a function of the center-of-mass momentum pjc+ of the K + meson at the 
Pk- values of 1.4 GeV/c and 1.8 GeV/c. (b) Relative contributions of phase-space (dashed line) and the 
production amplitude (solid line) parts to the cross section of the same reactions as in (a) at pk- =1-4 GeV/c 
for the same K + angle and the H mass. The curves are normalized to the same maxima. 



In Fig. 4(a), we show the double differential cross section dcr/dQ.K+dpK+ for the same reaction 
as in Fig. 3 at the beam momenta of 1.4 GeV/c and 1.8 GeV/c, corresponding to the CM angle 
6 K + = 0°, as a function of the CM momentum of K + (p K ¥)- The //-dibaryon mass was taken to 
be 2.224 GeV which is in the region where it is more likely to be located. It is seen that the cross 
sections are peaked very close to the kinematically allowed maximum regardless of the beam 
momentum. Since we have used the correct three-body phase-space factor, our K + spectra are 
broader than those given in Refs. Q20G and [21]. This is evident from Fig. 4(b) where we show the 
contributions of the phase-space only and the production amplitude only components to the K + 
spectrum at the beam momentum of 1.4 GeV/c and for the same value of m H as in Fig. 4(a). The 
curves are normalized to the same maxima so that their respective widths can be easily compared. 
We see that while the production amplitude part has a relatively narrower width and it peaks right 
at the maximum possible p K f, the three-body phase-space part is broad and has a maximum at 
somewhat lower p™. This results in a broader K + spectrum with a peak at somewhat below the 
maximum value of the K + CM momentum. It is likely that the use of the approximate phase-space 



factor of Refs. U 
authors. 



2 111 is the reason for the narrower widths of the K + spectrum obtained by those 



10 



In summary, we have studied the production of the stable six-quark //-dibaryon, which has 
spin-parity + , isospin and strangeness -2 via the (K~,K + ) reaction on a 12 C target within an 
effective Lagrangian model. The model assumes this reaction to proceed in two steps. In the first 
step, 3" hyperon and a K + meson are produced in the initial collision of the K~ meson with a 
target proton. In the second step, the E~ hyperon fuses with another target proton to produce the 
//-dibaryon. Although, similar models have been used earlier for such investigations, we have 
made several improvements to them. We employ a proper three-body phase-space to describe the 
A(K ,K + )HB reaction. The 3" + K + production amplitude has been calculated by excitation, 
propagation and decay of A and S hyperon resonance intermediate states in the initial collision 
of the K~ meson with a target proton. The vertex parameters (the coupling constants, and the 
form factors) at the resonance vertices have been taken to be the same as those fixed earlier by 
describing both the total and the differential cross sections of the elementary p(K + , K~)K~ reac- 
tion within a similar model. The same parameters were also used previously in the calculations 
of the 3 hypernuclei. The bound proton spinors have been obtained by solving the Dirac equation 
with vector and scalar potential fields having Woods-Saxon shapes. Their depths are fitted to the 
binding energy of the respective state. This is in contrast to the previous studies where S~ + K + 
amplitudes were obtained by parameterization of the not so well known experimental differential 
cross sections of the p(K + ,K~y£r reaction at 0°. Unlike the previous studies where the numeri- 
cal calculations were limited to a light 3 He target, we have applied our model to compute cross 
sections for a 12 C target. This makes it possible to make comparison, albeit indirectly, with the 
currently available limited data and to make predictions for future measurements. 

We find that the H dibaryon production cross sections in the 12 C(K~ , K + )H w Be reaction are 
more than an order of magnitude larger than those calculated earlier on a 3 He target. The mag- 
nitude of this cross-section for a H mass very close to the AA threshold, is comparable to the 
H production cross section determined in a study of the n C(K~ , K + )AAX reaction just above 
this threshold with a 90% confidence level. We notice that the differential cross-section for the 
12 C(K~, K + )H 10 Be reaction for observing K + at zero degrees peaks around the beam momentum 
of 1.4 GeV/c, which is above the production threshold of this reaction by almost the same amount 
as the position of the maximum is in the zero degree differential cross section in the elementary 
p(K~ , K + )"Br reaction above the corresponding threshold. 

The spectrum of the K + meson has a peak very close to the kinematically allowed maximum 
K + momenta and its width is broader than that observed in previous studies of this reaction on 
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a 3 He target. This is the consequence of using correct three-body phase-space and kinematics. 
The broad width of this spectrum might ease the experimental constraints on their measurements, 
with a caveat that it might become comparable with those of some background processes like 



2lh . However, the 



the momentum spectrum of K + recoiling against the continuum E~p pair 11201 
relatively large magnitude and shifted peak position of the cross sections for the latter process 
should make it possible to separate them from the H production events. 

This work was supported by the Dutch Research Foundation (NWO), and the University of 
Adelaide and the Australian Research Council through grant FL0992247(AWT). 
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